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ABSTRACT 

We have measured the mean pecuhar motions of 103 early- type galaxies in 7 clusters in 
the Perseus-Pisces (PP) ridge or PP background, and a further 249 such galaxies in 9 
calibrating clusters from the literature, using the inverse Fundamental Plane relation. 
This relation is found to have a distance error of 20% per galaxy. None of the 6 
clusters in the PP ridge has a significant motion with respect to the cosmic microwave 
background (CMB) frame, but the PP background cluster J8 shows marginal evidence 
of 'backside infall' into the PP supercluster. We find that the full 16 cluster sample has 
a mean CMB-frame bulk motion of 420±280 km s^^ towards I = 262°, b = -25°. This 
result is consistent both with no bulk motion in the CMB frame and with the ^ 350 
km bulk motion found by Courteau, Faber, Dressier and Willick. It is inconsistent 
at the 98% confidence level with the ~ 700 km s""'^ bulk flow found by Lauer & 
Postman. The PP ridge clusters are found to have a small and statistically insignificant 
mean radial motion with respect to the CMB frame: —60 ± 220 km s^^. Our error 
analysis fully accounts for the uncertainties in the mean Hubble flow, as well as the 
errors due to the merging of different spectroscopic datasets. A comparison between 
our cluster peculiar velocities and the predicted peculiar velocities from the IRAS 1.2 
Jy density field, smoothed on a 500 km s~^ scale, yields /9/ = Q^-^/bj — 0.95 ± 0.48, 
consistent with previous results. We find agreement between our peculiar motions and 
published TuUy-Fisher results for the same clusters. The disagreement between the 
11 clusters common to our sample and that of Lauer & Postman, based on brightest 
cluster galaxies (BCGs), is statistically significant at the ^ 99.7% confidence level 
indicating that the errors of one or both of these data sets are underestimated. When 
the BCG distances corrected for the X-ray luminosity of the host cluster are used, the 
disagreement is reduced to the '^94% confidence level. 

Key words: galaxies: distances and redshifts — galaxies: elliptical and lenticular, 
cD — galaxies: clusters: general — cosmology: observations — large-scale structure of 
Universe 



1 INTRODUCTION 

The two dominant concentrations of galaxies in the nearby 
Universe are the Great Attractor (hereafter GA) superclus- 
ter complex and the Perseus-Pisces (hereafter PP) super- 
cluster. The GA appears to be extended along the line of 
sight at I ~ 310°, b ~ 20°, with its densest region at and be- 
hind the Centaurus cluster in the distance range ~ 3000 - 
4500 km s~^ (Hudson 1993). On the opposite side of the sky, 
the main ridge of PP is extended perpendicular to the line 
of sight with 4500 ^ czq ^ 5500 km s~^ (Giovanelli and 
Haynes 1985; Wegner, Haynes and Giovanelli 1993). The 
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only way to measure the distribution of mass on these large 
scales is by measuring the peculiar velocities of galaxies and 
of clusters of galaxies. 

Studies of the peculiar flow field around these domi- 
nant concentrations can be traced from Shaya (1984) and 
Tammann & Sandage (1985) who found that, after having 
accounted for Virgo infall, nearby galaxies showed a resid- 
ual motion with respect to the frame of the Cosmic Mi- 
crowave Background (hereafter CMB) in the direction of 
Hydra-Centaurus. A tidal shear of the flow field in this di- 
rection was also noted by Lilje, Yahil & Jones (1986). Strong 
infall into a Great Attractor behind the Centaurus cluster 
at cz ~ 3000 km s"'^ was first claimed by Lynden-Bell et 
al. (1988). While it is now evident that there is a strong 
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streaming motion of galaxies in the direction of Centaurus, 
it remains unclear whether this motion is generated locally, 
e.g. by Centaurus or the GA, or whether it is due to more 
distant sources. For example, the claim by Mathewson, Ford 
and Buchhorn (1992) that 'backside infall' is not observed 
behind the Great Attractor, argues for coherent streaming 
motions generated at very large distances. 

In practice, the flow field is likely to be complex, be- 
ing dominated by attractors and voids on a range of scales. 
Nevertheless, to flrst order, the mean motion of the PP su- 
percluster allows a test of whether the flow is due to local 
or distant sources since, in combination with peculiar veloc- 
ity data near the GA, both the bulk motion and the shear 
of the peculiar velocity field can be measured. If the large 
bulk motion of nearby galaxies is due to local sources, then 
the peculiar velocity of PP will be small: ~ —100 km s~^. 
Alternatively, if the source of the motion is more distant, 
then PP should participate in a large-scale coherent bulk 
flow and have a peculiar velocity of ~ —400 km s~^ with 
respect to the CMB. 

Previous work on the motion of PP has been based 
mainly on the spiral galaxy TuUy-Fisher (TF) relation, the 
PP region not being well-sampled in the elliptical galaxy sur- 
vey of Faber et al. (1989). Willick (1990, 1991), reported a 
mean radial motion of —440 ±50 km s~^ for PP field spirals 
with redshifts 3800 < cz < 6000 km s~^ based on r-band TF 
data. The small quoted error on this result includes only the 
random error from the TF scatter. The result is subject also 
to a systematic calibration error of 2%, or 100 km s'^ at PP 
(Willick 1991). Courteau et al. (1993) added further r-band 
TF data and found that the PP region participated in a uni- 
form bulk flow of 360 ± 40 km s'^ towards I = 294° ,6 = 0°. 
Using I-band TF peculiar velocities of clusters in PP, Han 
& Mould (1992, hereafter HM) claimed a motion of -400 
km s~^. In contrast, the recent TF field survey of Giovanelli 
et al. (1996; da Costa et al. 1996) indicates no net motion 
of PP. 

The PP supercluster is at least as overdense in galaxies 
if not more so than the GA (Saunders et al. 1991; Hudson 
1993; Strauss & Willick 1995). If PP is as massive as implied 
by its overdensity in galaxies, then strong infall of galaxies 
into the PP supercluster is predicted. The infall signature 
is difficult to detect unambiguously using field samples be- 
cause individual galaxies scatter out of the filament, creating 
a spurious infall signature. This inhomogeneous Malmquist 
bias (Hudson 1994a; Dekel 1994) is particularly severe in 
the PP region due to the sharpness of the density con- 
trast between the supercluster filament and the foreground 
and background voids. While cluster distances are far less 
severely affected by Malmquist biases, the clusters typically 
reside along the ridge line of the PP filament itself, so that 
they are not efficient tracers of the infall pattern. 

In this paper, we use the Fundamental Plane (hereafter 
FP) distance indicator to measure the peculiar velocities of 
early-type galaxies in clusters in and behind the PP fila- 
ment. The use of a cluster sample minimises the Malmquist 
bias problems which plague field samples (Strauss & Willick 
1995). Furthermore, since the cores of clusters are rich in 
early-type galaxies, cluster-elliptical samples are less prone 
to contamination than spiral samples. 

The outline of this paper is as follows. In Section ^ 
we discuss the cluster and galaxy samples. In Section 0, we 



derive the FP relation for our combined data sample. In 
Section ^, we present the distances and peculiar velocities of 
individual clusters, and in Section |^ we model the fiow field 
and measure the bulk motion of the sample and the mean 
flow of the clusters in the PP filament. Section ^ discusses 
potential systematic effects. In Section |^ we compare our 
results to previous work. Finally, in Section ^ the paper is 
summarized. 

Note that distances - either the estimated distance d or 
the true position r - are quoted in units of km s^^. Radial 
peculiar velocities, u = v ■ r = cz — d, are with respect to 
the CMB frame unless otherwise noted. 



2 DATA 

The principal goal of this work is to determine the motion of 
the PP supercluster and compare this motion with the pre- 
dictions of flow models. In order to derive distances to the 
PP clusters we must adopt a distance indicator zero-point. 
The approach often used is to adopt one cluster, usually 
Coma, as a calibrator which is assumed to be at rest with 
respect to the CMB. However, individual clusters will have 
peculiar velocities with respect to the CMB of order ~ 300 
km s~^ (Grammann et al. 1995). For the case of Coma, this 
would translate to a zero-point uncertainty of 4%, and hence 
to an uncertainty of ~ 200 km s~^ in any derived motion of 
the PP supercluster. In order to reduce this uncertainty, we 
use here a set of clusters as calibrators and simultaneously 
fit for the distance indicator zero-point and a fiow model. 
To break the degeneracy between the zero-point and a bulk 
motion, we require a sample of clusters with good sky cov- 
erage. 

In Smith et al. (1997, hereafter PPI), we reported 
the spectroscopy and R-band photometry for early-type 
galaxies in seven clusters in the PP region. This includes 
the six clusters (7S21, Pisces, HMS0122, A262, A347 and 
Perseus [=A426]) in the PP ridge (4500 <, czq <, 5500 
km s~^) and the background cluster J8 at czq ~ 10000 
km s~^. Our 'calibration' clusters are from the published 
studies of Lucey and collaborators (Lucey & Carter 1988, 
hereafter LC88; Lucey et al. 1997, hereafter LGSC) and 
J0rgensen, Franx & Kjargaard (1995a, 1995b, 1996, here- 
after JFK95a, JFK95b, JFK96). We exclude the poor 
clusters from these studies, i.e. Pavo II, S639, S753, 
Doradus(=Grml3) and Grml5. The Centaurus cluster is 
also excluded. Centaurus has a complicated internal struc- 
ture, i.e. the Cen30 and Cen45 components, and is lo- 
cated in the GA foreground, so that it is not likely to 
be a good calibrator. This culling restricts our calibration 
sample to nine clusters, viz. A194, A539, Hydra(=A1060), 
Coma(=A1656), A2199, A2634, A3381, A3574(=K27) and 
DC2345-28(:=A4038=K44). Details of aU 16 clusters con- 
sidered here, including the data sources used, are given in 
Table |l| Note that in our analysis, we solve for the distances 
of all 16 clusters simultaneously and hence the term 'cali- 
bration sample' is somewhat artificial. 

For the Coma, A2199 and A2634 clusters, we use the 
spectroscopic and V-band photometric data as tabulated in 
LGSC. The spectroscopic data is from several systems, in- 
cluding a few measurements from our PP runs, but excludes 
the erroneous 'FLEX' velocity dispersion data of Lucey et al. 
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Table 1. Cluster Sample 
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(1991). The data for the remaining 6 rich clusters are from 
LC88 and JFK. For these six clusters we take the spectro- 
scopic measurements of LC88 and JFK95b, and the photom- 
etry from JFK95a's r-band measurements. Our full sample 
consists of 352 galaxies in 16 clusters. 

All spectroscopic data have been corrected for aperture 
effects following JFK95b and brought onto a common sys- 
tem using the derived offsets and weightings given in PPl. 
The effective surface brightnesses, {p)o, are k-corrected, cor- 
rected for + surface brightness dimming and for Galac- 
tic extinction. For the latter, we use the Burstein & Heiles 
(1982) values of E{B-V) multiplied by the factors 2.35, 3.02 
and 2.50 to obtain extinctions in the R, V and Gunn r bands, 
respectively. Surface brightnesses are transformed from the 
V and Gunn-r bands to the Kron-Cousins R band by sub- 
tracting 0.57 magnitudes and 0.33 magnitudes respectively 
(see PPI). Uncertainties in these transformations introduce 
only 0.5% distance errors and so are neglected. 

From all sources, we exclude galaxies with morpholog- 
ical type so/a and later or which are classified as morpho- 
logical rejects (type R) in PPI. The following background 
galaxies were also excluded: A01094 in A262 {cz = 14620 
km s"^), ZH56 in A194 [cz = 8269 km s~^) and RMH30 in 
Hydra {cz = 10672 km s~^). Finally, the galaxies RB28 in 
A2199, B03C in A347 and R338 in Hydra are outliers (at 
the > 3.5a level) in the FP relation found below. They are 
excluded from the analysis in this paper. 

The recession velocity adopted for each cluster is the 
mean redshift of the observed sample of galaxies in the clus- 
ter. In Table |l|, the parameter Ecz — CTci/xA^ gives the error 
in the mean sample redshift, where cjci is set to 500 km s~^ 
for clusters with n < 10 members, and is set to the mea- 
sured cluster dispersion or 500 km s~^, whichever is larger, 
for clusters with n > 10. We find that our mean sample 
redshifts and published mean redshifts for the same clusters 
agree to within ^ e^z- 

3 THE FUNDAMENTAL PLANE 

The methodology adopted in this paper is to fit cluster FP 
relations a priori without the use of redshift information. 



then to convert the cluster zero-points into cluster distances 
and finally to fit flow models a posteriori. Following Strauss 
& Willick (1995), we refer to this procedure as 'Method F. 
In this section we describe the FP fits. The Malmquist bias 
corrections and cosmological corrections needed to trans- 
form the FP zero-points into cluster distances are deferred 
to Section n. 



3.1 Method of fit 

The Fundamental Plane can be described by either a 'for- 
ward' or an 'inverse' fit, depending on whether the slopes are 
obtained by regressing on the distance dependent parameter 
— the logarithm of the effective radius, log Rc — or by re- 
gressing on the distance independent parameter, log a. Note 
that, in the inevitable presence of selection effects and scat- 
ter, the slope of the inverse relation is not simply the inverse 
of the forward relation slope. 

Galaxies in cluster samples are typically selected ac- 
cording to their morphological properties and either sur- 
face brightness, magnitude, diameter or some combination 
of these. The inverse relations have the advantage that they 
are unbiased by selection on the dependent variables. For 
the FP, this means that any selection in i?c, {/i)e or any 
function of these quantities (e.g. total magnitude) will not 
bias the inverse FP fit. 

For the PP clusters, galaxies were selected from pho- 
tographic plates, and are limited by APM or GSC photo- 
graphic magnitudes. From the candidate lists, spectra were 
obtained typically from brightest galaxy on down. There 
are, however, cases were the photographic magnitude was 
re-estimated by eye (e.g. due to contamination by a nearby 
star or galaxy) or where fainter galaxies were preferentially 
observed because they could be placed on the same spec- 
trograph slit as a brighter galaxy. Therefore, the selection 
criteria are likely to be somewhat fuzzy and are difficult 
to quantify precisely. As a result it is difficult to apply the 
forward-relation bias-correction scheme described by Willick 
(1994). Because the inverse distance indicator relations are 
insensitive to photometric selection, we shall consider only 
the inverse FP distance indicator in this paper. Of course 
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Table 2. Fit parameters and scatters for the distance indicator 
relations 



a 
P 

A, 



FP 

1.383 ± 0.040 
0.326 ± 0.011 

0.062 

0.198 



1.419 ±0.044 

0.065 
0.213 



the inverse relation is sensitive to explicit selection on ve- 
locity dispersion, but such selection is not present in our 
sample because we have not thrown out any galaxies a pos- 
teriori based on their velocity dispersions. There is one ex- 
ception to this rule: the Coma sample of LGSC explicitly 
excludes all galaxies with o < 100 km s~^. For this cluster 
we make a bias correction for the a cut which is exactly anal- 
ogous to the bias correction prescription of Willick (1994) 
for the calibration of the forward TuUy-Fisher relation from 
a magnitude-limited sample of cluster galaxies. In practice, 
this correction has only a very small effect due to the large 
range of o in the Coma sample. 

The inverse FP relation used here is: 

logo- = i log 7?o - -{m)c - -7ci (1) 
a a a 

We minimise the log a residuals over the galaxies in all 

clusters simultaneously, assuming the same slopes (a, /3)Q 

for all clusters but allowing individual zero-points (7ci) for 

each cluster to vary independently. In order to determine the 

errors on the zero-points, the slope of the FP relation are 

held fixed. Note that our slopes, derived from an inverse fit, 

will not be the same as those obtained from the "orthogonal 

fit" of JFK96. 



3.2 Results of FP fits 

The FP relation for the 16 cluster sample is shown in Fig. |l| 
Here, the mean CMB redshifts of the clusters have been used 
to shift all galaxies to the distance of the Coma cluster. 

Table gives the slopes, scatter in velocity disper- 
sion, Act, and the fractional inverse distance error, Ainv = 
a In(lO) Act, for the FP relation. Note that the distance er- 
ror per galaxy is larger than that found from the orthogonal 
fits of JFK96, through the influence of the larger slope, a. 
For comparison, we give also the slope and scatter of the 
inverse Dn — cr relation, for the same data. This distance in- 
dicator has slightly larger scatter than the FP, and is not 
used in the analysis to follow. 

Figures show the data for each cluster in separate 
panels. The solid line shows the distance indicator relation 
using the slope derived from all galaxies. The dotted line 
shows the median fit to the residuals. The difference be- 
tween the mean and median zero-point is always Si la, with 
the exception of A262 for which they differ at the 1.6a level. 
The scatter around the mean (solid) line is indicated in the 
lower right-hand corner of each panel. The individual clus- 
ter scatters are all consistent with the global scatter, with 
the exception of Coma, the scatter for which is marginally 



t Our definition of /3 differs by a factor -2.5 from that of JFK96. 
Here f} is the coefficient of (/i)e whereas for JFK96 it is the coef- 
ficient of log (7)g 



smaller: 0.046 versus the global scatter 0.062. This differ- 
ence is only marginally significant (at the ~ 95% confidence 
level). Note that since Coma is traditionally the standard 
cluster for FP/Dn — o studies, its galaxies have been ob- 
served many times. Thus their mean velocity dispersions will 
have smaller random errors than galaxies in less frequently 
observed clusters. The dashed line in each panel shows the 
individual fit with free slope to the cluster in question. The 
slopes (relative to the slope of best fit) are given in the top 
left hand corner. None of these slopes are significantly dif- 
ferent from unity, except for Pisces, for which the difference 
is marginal (~ 2cr). 



4 DISTANCES AND PECULIAR VELOCITIES 
4.1 Bias corrections 

We argued in Section ^ that the inverse fits were insensitive 
to selection on R^, {/i)o or functions thereof. As discussed 
above, we perform in this paper a Method I peculiar veloc- 
ity analysis. Whenever a Method I procedure is used, the 
distances will be affected by Malmquist bias, which depends 
on the density field of the peculiar velocity tracer popula- 
tion, before selection. For a Method I analysis using forward 
distance indicator relations, this would be the only source 
of bias. In contrast, a Method I analysis using inverse re- 
lations is affected, in general, both by Malmquist bias and 
by 'object selection bias'. For samples of field galaxies, the 
object selection bias is strong, particularly at the distance 
where the selection function peaks and beyond. For samples 
consisting of rich clusters, as is the case here, the probability 
that an entire cluster will drop out is a very weak function 
of distance and hence the object selection bias is negligi- 
ble, leaving only the Malmquist bias. For further discussion 
of these issues the reader is referred to Strauss & Willick 
(1995). 

One of the advantages of a cluster sample is that the 
corrections for Malmquist bias are smaller than for field 
galaxies. The fractional distance error for each cluster is 
Ainv/\A'') where n is the number of observed galaxies in 
the cluster. The homogeneous Malmquist correction, which 
assumes that the underlying density field of clusters is uni- 
form, is a multiplicative factor exp(— 3.5Afjj^/n). For our 
sample, the largest homogeneous Malmquist correction is 
2% (or ~ 100 km s~^ at the distance of the PP ridge) for 
the cluster 7S21 which has 7 observed members. The correc- 
tion for the Perseus and Pisces clusters, which dominate the 
fiow by virtue of their small random errors, is only ~ 0.5% 
or ~ 25 km s~^. The correction due to the fact that the 
underlying cluster density field is inhomogeneous (Hudson 
1994a; Dekel 1994; Strauss & Willick 1995) is expected to be 
smaller than the homogeneous Malmquist correction, which 
itself is quite small, as demonstrated. We therefore neglect 
the inhomogeneous correction in this paper. 



4.2 Distance calibration 

In order to obtain cluster distances, we shall initially fix the 
distance and peculiar velocity of Coma to be 7200 km s~^ 
and zero, respectively, and calculate cluster distances rela- 
tive to Coma. In the fiow model fits to follow, we will allow 
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Figure 1. The FP relation for all galaxies. All clusters are shifted by their relative mean CMB redshifts, so that they are placed at 
the distance of Coma. Early type galaxies (E, E/SO, D or cD) are indicated by circles, later types are indicated by triangles. The 103 
galaxies in the 7 PP clusters are represented by filled symbols. The remaining 249 galaxies in the 9 calibrating clusters are shown by the 
open symbols. The rms scatter in logc, Ao-, for the full sample, and the E and SO s ubsamples is given in the lower right hand corner, 
with the offset between E and SO galaxies in the inverse FP relation. See Section 



5 for further discussion. 



this zero-point to be rescaled by a free parameter. We shall 
see below that this free parameter is close to unity, thus 
justifying our calibration a posteriori. 

The difference between the zero-point of a cluster and 
that of Coma is the logarithm of relative apparent angular 
diameters of galaxies (at the same a). We iteratively solve 
the angular diameter distance equation with go = 0.5 in 
order to convert the zero-point offsets to distances, i.e. the 



cz that the cluster would have if it were following the Hubble 
flow. For qo = 0, the distance of clusters in the PP ridge 
changes by only 10 km s~^. 

4.3 Results 

The Malmquist-corrected FP distances, diviCi and radial pe- 
culiar velocities, uqmb, are given in Table BI The distance 
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Figure 2. FP data and fits for the PP clusters 7S21, Pisces, HMS0122, A262, A347, Perseus and J8. Early type galaxies (E, E/SO, D 
or cD) are indicated by circles, later types are indicated by triangles. The solid line shows the global inverse FP, found by minimising 
logo" residuals simultaneously over the whole cluster sample with the same slope but varying zero-points for each cluster. Each cluster's 
measured scatter in logo" around this global fit is given in the lower right-hand corner. Galaxies which deviate from the global fit by 
more than 2.5 times the global scatter are labelled. The dotted line shows the median of the residuals with the slope fixed from the whole 
cluster sample. The dashed line shows the best slope and zero-point fit to the individual cluster. The individual cluster slope relative to 
the global slope is given in the upper left-hand corner. 



error, e^, is due to the error in fitting the zero-point. The 
pecuhar velocity error, £„, is the distance error added in 
quadrature -with the cluster mean redshift error, Ecz. 

Figures ^ and ^ show the CMB redshifts and the FP- 
inferred distances of individual galaxies in each cluster. Note 
that there is no tendency for galaxies to lie along the Hub- 



ble line, in contrast -with some spiral-cluster samples (Willick 
et al. 1995). By observing early- type galaxies, "we have effi- 
ciently selected galaxies in the virialized cluster cores. 

The Hubble diagram for all 16 clusters in shown in 
Fig. 0. Note that most clusters have a peculiar velocity in 
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Figure 3. FP data and fits for the clusters A2199, A2634, Coraa, A194, A539, A3381, A3574, DC2345-28 and Hydra. Symbols and 



curves are as in Fig. ^ 



the CMB frame which is not different from zero by > 2a. 
The exceptions are A2634 and A194. 

Figures ^ and ^ show the flow field of the PP clusters 
studied here. The large circles show the FP distance to each 
cluster, while the tip of each arrow lies at its mean redshift 
in the CMB frame. The length of the vector thus gives the 
peculiar velocity. The horizontal lines bracketing each clus- 
ter position indicate the distance errors. The small points 
are plotted at the (redshift-space) positions of galaxies from 
the CfA ZCAT compilation (Huchra et al. 1993), and are 
shown to illustrate the PP structure as traced by galaxies. 



It is clear from these figures that the mean CMB-frame mo- 
tion of the PP ridge is small, and that the motion of PP 
ridge clusters with respect to the mean motion of the ridge 
is quiet. Note also that the cluster J8 shows some evidence 
of backside infall into the PP supercluster. 



5 THE FLOW FIELD 
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Figure 4. Distances and redshifts of individual galaxies in 7S21, Pisces, HMS0122, A262, A347, Perseus and J8. Horizontal and vertical 
solid (dotted) lines give the mean (3it range) of redshift and log(distance), respectively. The short dashed lines give the respective 
medians. The Hubble flow is shown the diagonal long-dashed line. The offset of the cross from the Hubble line gives the cluster peculiar 
velocity. 



5.1 Flow models 

In order to model the peculiar velocities, we will consider 
flow models composed of a combination of pure Hubble flow, 
a uniform bulk flow and/or the peculiar velocity predicted 
from the IRAS 1.2Jy redshift survey (Fisher et al. 1995a). 

In order to allow for a possible 'thermal' component 
in the peculiar velocity of Coma, our adopted calibrating 



cluster, we allow the zero-point of the distance scale to float 
by including a monopole (or Hubble) term 

VH{r) = {I + An)r . (2) 

The bulk flow is modelled as a uniform dipole motion, 
independent of distance: V . 

If light traces mass, then the most reasonable flow mod- 
els are those in which the galaxy density fleld is used to pre- 
dict peculiar velocities. There is usually one free parameter 
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Figure 5. As Fig. | for tlie clusters A2199, A2634, Coma, A194, A539, A3381, A3574, DC2345-28 and Hydra. 



for these models, /3 = /(n)/b ~ ^^'^ jh. This allows us to 
fit a degenerate combination of il, which affects the pecu- 
liar velocity field for a given mass fluctuation field, and the 
biasing parameter 6 = (5gai/5mass which relates the observed 
density fluctuation field of galaxies to that of the mass. Note 
that different tracers of the galaxy density field (e.g. IRAS 
galaxies versus optically-selected galaxies) will have different 
bias parameters and hence different values of /9. 

In this paper, we use predicted peculiar velocities de- 
rived from the IRAS 1.2 Jy redshift survey density field, 
which was kindly provided to us by M. Strauss. The IRAS 
density field extends to 12000 km s~^, which allows us to 



make predictions for all of the clusters in our sampler! Rich 
clusters of galaxies have collapsed from scales ~ 8/i~ Mpc, 
so we expect that their peculiar velocities are accurately de- 
scribed by the predictions of linear theory: 



-d r . 



(3) 



where &i is the IRAS density fluctuation field, smoothed 



t All-sky density fields of optical galaxies (e.g. Hudson 1993; San- 
tiago et al. 1995 ) are not sufficiently deep for our purposes. 
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Figure 7. Flow field for PP clusters with the open circle indicating the FP distance of the cluster and the straight lines showing the 
1(7 distance error. The tips of the arrow show the mean cz in the CMB frame. Also displayed are the CMB redshift space positions of 
galaxies taken from the CfA ZCAT compilation. Note that the 'cluster' of points just west of J8 is a thin wall which runs north-south 
and is seen here in projection. 



with a 500 km radius Gaussian filter and processed 
with a Weiner filter in order to reduce the effects of shot 
noise. The smoothing scale is approximately equivalent to 
an 8/i~^Mpc top-hat, and so is a good match to the scale 
from which the clusters collapsed. Note that linear theory 
has been used to transform iteratively the IRAS galaxy po- 
sitions from redshift space to real space. The /?/ adopted 
for this iterative scheme is 1.0, which is not significantly 
different from the value obtained from fits to our peculiar 



velocity data below. Note that the effect of the /9 used in the 
redshift-to-real space iteration procedure on the derived /3 

is typically very small (Hudson 1994b; Hudson et al. 1995). 

Errors in the IRAS predicted peculiar velocities can be 
estimated from the results of Fisher et al. (1995b). Those 
authors generated mock IRAS 1.2 Jy surveys from N-body 
simulations. Although they used a different redshift-to-real 
space reconstruction scheme (based on a spherical harmonic 
decomposition) , in common with the case here they also used 
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Figure 8. As in Fig. M, but only out to a distance of 6000 km s , to show more detail in the PP ridge. 



the Weiner filter to suppress shot noise. We therefore expect 
their IRAS predictions to have similar error properties to 
those used here. Fisher et al. estimated that the errors in 
the IRAS-predicted radial peculiar velocities of individual 
galaxies were ~ 200 km s"^ within 6000 km s"\ growing to 
approximately 250 km s~^at 11000 km s~^. We expect that 
the errors for predicted peculiar velocities of clusters should 
be somewhat smaller than these values. 

In general, the predicted peculiar velocity is then 

Vp{r) = VB_{r) + V + vi{r) . (4) 

Note that for those cases in which we use the IRAS predic- 



tions as part of the flow model, V should be interpreted not 
as the mean bulk flow of the sample, but as a residual bulk 
flow due to mass fluctuations beyond the limit of the IRAS 
density fleld. 

Given the model predictions at the estimated distances 
of each cluster, we minimise, with respect to the free param- 
eters Ah, V and Pi, the statistic 



(5) 



where eth allows for an additional 'thermal' peculiar veloc- 
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Table 3. FP Cluster Distances and Peculiar Velocities 



Cluster 












km s~ "'^ 


km 


km s^ ^ 


km s "'^ 


7S21 


5517 


5448±409 


69±450 


-201 


Pisces 


4714 


4583±182 


131±208 


108 


HMS0122 


4636 


4680±310 


-44±352 


159 


A262 


4528 


4787±300 


-259±340 


226 


A347 


5312 


5589±392 


-277±430 


-210 


Pgj-geus 


5040 


5176±185 


-136±307 


180 


J8 


9425 


10457±576 


-1032±602 


-706 


A2199 


8947 


9289±307 


-342±325 


-29 


A2634 


9158 


10118±339 


-960±364 
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7200 


7200±170 


0±204 
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A194 
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295 
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8615 
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644 
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11471 
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A3574 
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Figure 6. Hubble diagram for all 16 clusters. 



ity of the cluster with respect to the flow model. We set 
Eth = 200 km s~^, so that the reduced values are ap- 
proximately unity. Note that we neglect errors in the IRAS 
predicted peculiar velocities, but these are expected to be 
< 200 km s-i. The error in the measured radial peculiar ve- 
locity is typically twice this value; we therefore expect that 
the derived /3/ will not be significantly biased if errors in the 
IRAS predictions are neglected. 

The formal random errors are obtained from the co- 
variance matrix in the usual way. This procedure, however, 
neglects an important source of error, namely the system- 
atic effects introduced by errors in matching the different 
velocity dispersion systems. In PPI, we determined the cor- 
rections which were applied to the raw log a measurements 
to bring them onto a common system, and described the 
bootstrap resampling method which was used to determine 
the uncertainty in these corrections. In order to quantify 



the effect that these uncertainties have on our results, we 
perform the following procedure. For each of 200 bootstrap 
samples, we obtain a realization of logo- offsets. These are 
then used to generate new merged spectroscopic data sets 
for the 352 galaxies in this paper. The new data are then 
passed through the flow analysis above. In this way, we de- 
termine the covariance matrix which is due solely to system 
matching errors. This covariance matrix is then added in 
quadrature to the random error covariance matrix to obtain 
the total error covariance matrix. The bulk motion errors 
quoted in the following subsection include these systematic 
contributions. 



5.2 Results 

We fit the bulk fiow and IRAS-predicted flow models to two 
samples of clusters: the full sample of 16 clusters, and the 
6 PP ridge clusters alone. Table ^ lists the values of /3, Ah 
and the amplitude and direction of bulk flow V for the full 
cluster sample. For the flrst two parameters, entries without 
a quoted error indicate that the value was held fixed in the 
fit. The of the fits and the number of degrees of freedom 
are given in the last column. 

When the covariance matrix of the bulk fiow compo- 
nents is diagonalized, we find that one direction has much 
smaller errors than the others. Hereafter this direction will 
be referred to as the DME for "direction of minimum error" . 
This effect occurs because the fits are dominated by the PP 
clusters, which are concentrated in a small patch on the sky. 
The rows of Table ^ give, for each solution in the correspond- 
ing row in Table ^, the directions of the three eigenvectors of 
the covariance matrix and the projections of the bulk mo- 
tion along the these directions. The last columns give the 
significance of bulk fiow compared with zero: x'^ (to be com- 
pared with three degrees of freedom) and the corresponding 
probability that the there is no bulk motion in the CMB 
frame. The first eigenvector, denoted by the subscript "1" is 
the DME. For the fuU cluster sample, the DME is / = 351°, 
h = 51°, which is within 20° of the of the bulk motion found 
by Lauer & Postman (1994, hereafter LP) and is close to 
the directions of the GA (Faber & Burstein 1988) and the 
Shapley Concentration (Raychaudhury 1989). Thus, despite 
the large errors in the other eigendirections, this sample can 
be used to test the predictions of a number of interesting 
fiow models. 

The first row of Table ^ indicates that the mean bulk 
motion of the entire sample with respect to the CMB frame 
is 420 ± 280 km s"^ towards / = 262.6°, b = -25.3°, which 
is not significantly different from zero. 

The fits with a free Hubble term. Ah, (Rows 1 and 3 in 
Table ^ indicate that the best fitting Ah is very small and 
consistent with zero. Thus, by chance, it would appear that 
Coma has a small peculiar velocity with respect to the best 
fitting flow model. We therefore flx Ah = for the fits to 
the PP ridge sample. 

For the PP ridge sample, the axis of the DME runs be- 
tween I ^ 314°, b ^ 26° and I ^ 134°, b ^ -26°, i.e. along 
the mean radial direction of the PP clusters, near A262. 
The total peculiar velocity error of 220 km s"^ in this di- 
rection can be broken down as 161 km s~^ (random error), 
125 km s~^ (system matching error) and 85 km s~^ (the er- 
ror arising from uncertainty in the distance scale zero-point 
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Table 4. Bulk flow fits to the sample of all 16 clusters 



Row 


13 
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1^1 


I 


b 


xVd-o.f. 
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0.00 
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0.001 ±0.017 


389 ± 294 


313.2 
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Table 5. Bulk flow eigenvalues and eigenvectors 
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Figure 9. FP radial peculiar velocities for clusters within 30° 
of the direction of minimum error (DME), I = 311°, b = 51°, as 
a function of their distance from the LG. The distances and pe- 
culiar velocity of clusters on the Perseus-Pisces side of the sky, 
i.e. opposite to the DME, have been multiplied by —1. Symbol 
area is proportional the effective weight of the cluster in the flt. 
The hatched regions shows the ±lcr range of peculiar velocities 
within the 30° semi-angle cone predicted by various flow models; 
the heavy solid line shows the mean predicted peculiar velocity 
within the cone. The flow models are labelled as followed: LP 
- the bulk motion of Lauer &: Postman (1994); GA - Faber & 
Burstein (1988). 



from column 2, row 1 of Table ^). Thus the system-matching 
error accounts for 18% of the total error. The errors on the 
components of the mean motion in the two other transverse 
directions are very large (^ 700 km s"^ and ^ 7000 km s"^), 
so the motion of the PP ridge is only well-determined in the 
radial direction. The mean radial peculiar motion of the PP 
ridge clusters is — 60±220 km s~^. The PP ridge is thus con- 
sistent with having no radial peculiar motion in the CMB 
frame. 

Fig. ^ shows the peculiar velocities projected along the 
DME as a function of their distance from the LG (with the 



PP direction taken to be negative and the GA direction 
positive). Only clusters within 30° of the DME are shown. 
The predictions of the GA, LP and IRAS flow models are 
shown schematically by the hatched regions, which indicate 
the ±1(7 range of predicted peculiar velocities within the 
cone. Note that the GA and IRAS predictions overlap in 
the region where the GA model was first defined, —2000 < 
d < 4000 km s"\ The IRAS predictions are more complex 
elsewhere, in particular infall around PP is expected. Note 
that the PP ridge clusters do not lie in the foreground or 
background infall regions, but right along the centre of mass 
of the supercluster where vi ~ 0. Finally, note that this 
diagram is intended to show the qualitative behaviour of the 
data and fiow models near the DME. It should not be used 
to compare the observed and predicted peculiar velocities of 
individual clusters since the predictions are only the mean 
within the cone and not the value at the position of the 
cluster itself. 

We now consider in more detail the predictions of the 
IRAS density field. Fig. ^ plots the observed peculiar veloc- 
ities of individual clusters against the IRAS predictions at 
their estimated positions, with /3/ = 1. The slope of the line 
of best fit thus gives the value of /3/: the result is 0.95 ±0.48. 
This is a marginally significant detection of /?/, but the er- 
rors are too large to place useful constraints on its value. In 
principle, however the comparison of cluster peculiar veloc- 
ity fields with the predictions from a deep all-sky redshift 
survey should lead to a determination of (3i which is more 
reliable than most other methods. The advantage of our ap- 
proach is that inhomogeneous Malmquist bias, smoothing 
effects and the effects of non-linear bias and gravity are all 
small for cluster peculiar velocities. 

For the fits with f3i 7^ 0, the quoted bulk motion is the 
residual CMB-frame motion due to sources not modelled 
by the IRAS density field, i.e. beyond 12000 km s~^ or due 
to shot noise errors within that distance (see Willick et al. 
1997b for a discussion of these errors). It is interesting to 
note that the prediction for the error-weighted mean mo- 
tion of the PP ridge clusters is 60/3/ km s~^ away from the 
LG. This has a direction opposite to the prediction of the 
GA model. Presumably the effects of other structures, such 
as the foreground void, are important for the dynamics of 
the PP ridge. The residual motions of both the full sam- 
ple and the PP ridge sample are small and consistent with 
zero. Finally, we note that the IRAS-predicted peculiar ve- 
locity for J8 is —706 km s~^ in good agreement with the 
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Figure 10. Observed and IRAS-predicted peculiar velocities of 
individual clusters. The predictions are scaled to /3j = 1 so the 
slope is a measure of /3j. Clusters with observed or predicted 
peculiar velocities larger than 500 km s~^ are indicated. The solid 
line shows the slope of best fit, and dashed lines give the Itr error 



observed value —1032 ± 602 km s~^. However, for A2634, 
the IRAS prediction is only —127 km s~^, which appears 
to be marginally inconsistent with the observed —960 ± 364 
km s~^. Note, however, that the shot-noise error in the IRAS 
prediction is ~ 250 km s~^ and the systematic uncertainty 
in the peculiar velocity of A2634 is ~ 230 km s~^. 

We have considered fits to the GA model of Faber & 
Burstein (1988). Our sample is not ideal, however, for con- 
straining the GA model, since it contains few clusters in the 
regions where the GA model predicts strong infall. In partic- 
ular the Cen30 and Cen45 clusters have been excluded from 
our sample. In performing fits of this type, we find poor 
agreement with the model. The disagreement does not arise 
from the small infall motion of the PP ridge, which is in 
rough agreement with the GA model, but is instead partly 
due to A3574 which has a peculiar velocity of 655 ± 369 
km s~^ near the centre of the GA, where its peculiar ve- 
locity should be ~ km s~^. This result should be taken 
with some caution since the peculiar velocity of A3574 may 
be biased by a low a galaxy which is an outlier in the FP 
relation (see Section 3.1 below). 

We have also considered models involving infall to the 
Shapley Concentration (Raychaudhury 1989). We find that 
the results are statistically indistinguishable from a pure 
bulk motion in the same direction, because we have no clus- 
ters close enough to the Shapley Concentration for its tidal 
shear to be important. 

Finally, we can use our results to constrain the thermal 
component of the clusters' peculiar velocity with respect to 
a given flow model. We do this by varying eth and requir- 
ing that the resulting (equation (^)) lie in an acceptable 
range. For the full cluster sample, we find that eth < 550 



km s ^ at the 95% confidence level, with the most likely 
value being ~ 250 km s-\ This agrees well with the typical 
values seen in N-body simulations (Gramann et al. 1995). 
The fiow field of the 6 clusters in the PP ridge seems to be 
particularly cold: the 95% upper limit on their rms peculiar 
velocity is ~ 150 km s~^. 



6 TESTS OF POTENTIAL SYSTEMATIC 
EFFECTS 

Systematic errors are often a major concern for distance 
measurements. In the analysis above, we have considered 
the uncertainty introduced by errors in matching the differ- 
ent velocity dispersion datasets. In this section we explore 
other effects that might lead to biases in our peculiar veloc- 
ity results. 



6.1 Effect of low-dispersion galaxies 

Galaxies with velocity dispersions below 100 km s~^ can 
have large systematic and random errors in the measured 
velocity dispersion (JFK95b). Our PP target sample con- 
tains only 5 galaxies with velocity dispersions smaller than 
100 km s~^, and the full sample of 16 clusters contains only 
16 such galaxies. We have tested the effect these galaxies 
might have on our solutions by explicitly excluding them, 
and correcting for the a cut in a manner exactly analo- 
gous to the bias correction prescription of Willick (1994) 
for the calibration of the forward TuUy-Fisher relation from 
a magnitude-limited sample of cluster galaxies. We find that 
none of our results (the slopes of the FP relation, the pecu- 
liar velocities of individual clusters or the global fiow model 
fits) change by an amount larger than the random error, 
with the exception of the peculiar velocity of A3574, which 
drops from 655 km s~^ to 66 km s~^ when the galaxy W69 
is excluded. 



6.2 Effects of galaxy morphology 

We have assigned galaxies to morphological categories us- 
ing their appearance on our CCD frames, and their radial 
profiles derived from the aperture photometry. The classifi- 
cations for PP galaxies is given in PPI. Note that the cat- 
egories are broad: E ; SO; Q (unclassified or unclassifiable; 
there are only two galaxies with this classification). 

To assess the validity of the equal treatment of E and 
SO galaxies in our analysis, we can compare the relative off- 
set of elliptical galaxies (E, E/SO, cD or D) and SOs (SO, 
SO/E) with respect to the mean inverse FP relation. Ellipti- 
cal galaxies have slightly higher log a at the same values of 
i?o-0.326{/i)o than SO types by 0.018±0.007. Using a sample 
which has a large number of galaxies in common with that 
considered here, JKF96 found a difference of 0.006±0.011 in 
the opposite sense. Our sample includes about twice as many 
galaxies as that considered by JFK and we have adopted an 
optimal correction scheme for combining the different ve- 
locity dispersion datasets. The scatter, Ao-, of the elliptical 
galaxies is 0.054 compared to 0.073 for the SO galaxies. This 
difference may be due to the fact that SO galaxies tend to 
have lower velocity dispersions which carry with them larger 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy clusters in the Perseus-Pisces region — //. The peculiar velocity field 15 



observational errors. We conclude that there is marginal evi- 
dence of a difference between FP relations of E and SO galax- 
ies. We have investigated the effects of this difference on our 
flow results and find that, because most clusters have simi- 
lar proportions of E and SO galaxies, none of our results are 
sensitive to this offset. 



6.3 Stellar population effects and the FP-Mg2 
relation 

The presence of a small percentage of stars with interme- 
diate ages would increase the scatter of the FP relation, 
and cause systematic biases in the derived distances. Since 
the distance- independent Mg2 index is also sensitive to stel- 
lar populations, the use of Mg2 as an additional parameter 
in the FP relation had been advocated (Guzman & Lucey 
1993). There are, however, several practical reasons why we 
prefer not to use Mg2 in our distance indicator. We shall 
see that for inverse fits, including Mg2 decreases the scatter 
in logo" but actually increases the random distance error, 
Ainv due to the change in slope a. Furthermore, by includ- 
ing Mg2 , the total system matching distance errors increase 
(by about a third) due to the addition of the Mg2 system 
matching errors and due to the velocity dispersions system 
matching errors, again because of the change in slope a. 
Finally, only 315 of the 352 galaxies in our sample have 
Mg2 measurements, and we are required to drop the cluster 
DC2345-28 which has only one Mg2 measurement amongst 
the spectroscopic systems considered in this paper. 

Nevertheless, in order to investigate the effects of stellar 
populations on our results, we have proceeded with fits to 
an FP-Mg2 relation of the form 

logcr = - log7?e - -{^i}c - -Mg2 - ^7,1 . (6) 
Q a a a 

We find the results a = 1.795 ±0.083, P = 0.331 ±0.016 and 
C = -2.38 ± 0.22 with = 0.054 and Ainv = 0.22. This 
highly significant detection of the Mg2 correlation arises as a 
result of the inverse nature of the fit. The log cr-Mg2 relation 
has a scatter of only 0.086 in log a, which is only slightly 
larger than the scatter in the inverse FP. The inverse fit 
thus gives considerable weight to Mg2 as a predictor of log a, 
in contrast to the forward FP. The distances of individual 
clusters typically change by ^ Icr, except for Coma which 
shows an offset at the 2.2a level. Because Coma was our 
nominal calibrator, all distances require rescaling by 5% with 
the result that Coma has a peculiar velocity of ±381 km s~^ 
and the remaining clusters follow the Hubble flow in the 
mean. The best fltting bulk flow of the full cluster sample 
is in fact somewhat smaller than that of the FP solution 
but well within the la errors. The motion of the PP ridge 
towards the LG increases slightly to 141 km s~^, which again 
differs by < la from the FP solution. 

We conclude that the residuals of individual galaxies 
from the inverse FP relation are weakly correlated both with 
galaxy morphology and with Mg2 index. However, when av- 
eraged over all galaxies in a given cluster, there is a negligible 
net effect of such correlations on derived cluster distances, 
and hence the flow results are robust against such effects. 



7 COMPARISON WITH OTHER WORK 

7.1 Bulk flow and the motion of Perseus— Pisces 

For the full 16 cluster sample, we find a bulk motion which 
is consistent with these clusters being at rest in the CMB 
frame. In order to compare our bulk motion, Vi, with co- 
variance matrix Ci with an independent sample with bulk 
motion V2 and covariance matrix C2, we calculate 

Xbuik = {Vi - ^2)^ (Ci ± C2)-' (Vi - V2) (7) 

and compare the result to a Xbuik distribution with 3 de- 
grees of freedom (following Hudson & Ebeling 1997, here- 
after HE)^ We find that the bulk motion of our sample 
is consistent with the 360 ± 40 km s~^ motion found by 
Courteau et al. (1993). However, the bulk motion found in 
this paper is inconsistent at the 98% confidence level with 
the result of LP, who used the photometry of brightest clus- 
ter galaxies (BCGs) as a distance indicator. Our result is 
also inconsistent at the 94% confidence level with bulk mo- 
tion obtained by HE, who applied a correction for the X-ray 
luminosity of the host cluster to the BCG distance indicator 
of LP. 

We find a negligible mean radial motion for the PP ridge 
(—60 ± 220 km s^^), which is, at face value, marginally in- 
consistent with the result of Willick (1990) who found a 
mean peculiar velocity of —441 ± 49 km s~^ (random error) 
from a TF study of PP field spirals in the redshift range 
3800 < cz < 6000 km s~^. Note, however, that the system- 
atic calibration error on this result is ~ 100 km s~^ (Willick 
1991). Furthermore, the two samples probe different parts 
of the PP supercluster: the PP ridge clusters studied here 
range from O'' ^ RA ^ 3*^ whereas the volume probed by 
Willick covers a wider range of redshifts and extends farther 
to the west, from 22'' ^ RA ^ 3*'. Willick (1991) noted that 
the western half of his sample had a greater infall toward the 
LG. It is not clear then that the two results are in significant 
disagreement. 

Our result is also marginally inconsistent with the result 
of HM who found a PP radial motion of -400 km s~^ from 
a TF cluster survey. However, it should be noted here that 
the HM and Willick (1990, 1991) results are not indepen- 
dently calibrated. Both studies use Aaronson et al. (1986) 
linewidths for the calibrating cluster sample, and linewidths 
from the work of Giovanelli, Haynes and collaborators (e.g. 
GiovaneUi & Haynes 1985; GiovaneUi et al. 1986) for the PP 
sample. HM also adopt the prescription of Willick (1991) to 
transform between the two systems. The transformation is 
based on a sample of 58 galaxies in common between the 
two linewidth sources. 



§ Note that when comparing two peculiar velocity samples with 
different sky coverage and effective depth, we do not expect the 
measured bulk flows, Vi and V2, to be identical even in the limit 
of no measurement errors due to the different window functions 
(Watkins & Feldman 1995). Therefore, the confidence with which 
we conclude that two samples are inconsistent will in general be 
slightly overestimated. 
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7.2 Individual cluster distances and peculiar 
velocities 

In order to understand the above-mentioned conflicts be- 
tween our results and those of Han & Mould, Lauer & Post- 
man and Hudson & Ebeling, we now turn to a comparison 
of distance and peculiar velocities of individual clusters. 

We have compared our cluster distances and peculiar 
velocities with the inverse TF cluster distances and peculiar 
velocities of HM, as rederived by Willick et al. (1997a) for 
the 7 clusters in common to the two samples. The results are 
shown in the upper panels of Fig. ^ and peculiar velocities 
are tabulated in Table ^ 

The two sets of cluster distances are marginally incon- 
sistent at the 90% confidence level. However, we find that 
there is good agreement between cluster peculiar velocities. 
Note that elliptical and spiral galaxies in the same nominal 
cluster do not necessarily sample the same kinematic ob- 
ject. In particular, spiral cluster samples typically extend to 
larger radii, and may be more prone to contamination from 
the field, or from infalling groups. 

For example, for the most discrepant cluster, A262, 
we find cz = 4528 km s~\ d = 4787 ± 300 km s-\ and 
hence u = —259 km s~^ whereas Willick et al. (1997a) 
find cz = 5057 km s"\ d = 5840 ± 322 km s"^ and hence 
u = —783 km s~^ Furthermore, the nominal cluster cen- 
tres differ by 1°. It is not clear, therefore, whether our A262 
cluster is the same object as that of HM. A similar situa- 
tion holds for the next most discrepant cluster, A2634, in 
which the treatment of the companion cluster A2666 affects 
the comparisons. Scodeggio et al. (1995) present a full dis- 
cussion of the contamination problems for this cluster. If we 
remove either one of these two clusters from both samples, 
the distances are compatible. 

In contrast to the TF results of Wilhck (1990,1991), 
Courteau et al. (1993) and HM, the recent I-band TuUy- 
Fisher cluster (Giovanelli et al. 1997a, b, hereafter SCI) and 
field (Giovanelli et al. 1996) data suggest no net motion 
for PP. The lower panels of Fig. ^ compare our distances 
with the cluster distances of Giovanelli et al. (1997a, b). We 
have computed the mean CMB redshifts of cluster mem- 
bers from Giovanelli et al. (1997a) using galaxies denoted 
by "c" in column 7 of their Table 2, corresponding to their 
"in" samples) and peculiar velocities determined from their 
incompleteness-corrected TF magnitude offsets from column 
5 of Table 3 of GiovaneUi et al. (1997b). Note that their NGC 
383 Group corresponds to Pisces and their NGC 507 group 
corresponds to HMS0122. While the two samples of cluster 
distances disagree at the 99% confidence level, this disagree- 
ment is mainly due to A2634 for which the mean redshifts 
differ by 650 km s~^. If this cluster is removed from the 
comparison, the distances are consistent. 

The data used to determine the peculiar velocities of 
A2199 and A2634 are the same as that used by LGSC. 
Whereas we derive values of — 342±325 km s~^ and — 960± 
364 km s~^ respectively, LGSC derived -160 ± 380 km s"^ 
and —670 ± 490 km s~^. These relatively small differences 
arise from the application of the FP distance indicator. In 
this paper we have used the 'inverse' relation whereas LGSC 
used the 'forward' relation with an allowance for the differ- 
ent selection functions in the clusters. 

Scodeggio, Giovanelli & Haynes (1996) have recently re- 



ported new FP results for Coma and A2634. Using Coma as 
the calibration cluster, they derive a distance for A2634 of 
9099 ±266 km s~^. While their photometric parameters are 
from a new set of I-band data, they used velocity dispersion 
data from several sources. For the Coma cluster these are 
mostly taken from the literature, while for A2634 they use 
mainly new velocity dispersion measurements from the Hale 
5m telescope. They have only a small overlap between their 
measurements and the literature values. Their quoted er- 
ror in the A2634 distance does not include the uncertainty 
in linking their velocity dispersion values onto the 'litera- 
ture' system. In PPI we showed that different spectroscopic 
datasets typically differ at the level of 0.01 dex. Using the 
new measurements from PPI and LGSC, we find that the 
Scodeggio et al. velocity dispersions are ~ 0.02 dex smaller 
than our standard system. While a more extensive analysis 
is required to place the Scodeggio et al. data accurately onto 
a 'standard' system, such a systematic offset would trans- 
late into increasing their derived distance for A2634 by ~500 
km s~^. After allowing for this correction the Scodeggio et 
al. distance for A2634 is not in conflict with the value re- 
ported here, i.e. ~ 9600 ± 350 km s"^ versus 10118 ± 339 
km s^^. 

Our non-detection of a bulk motion is in disagreement 
with the result of LP, who studied BCGs in 119 AbeU/ACO 
clusters within cz < 15000 km s~^. They used a distance 
indicator based on the photometry of BCGs and concluded 
that their sample had a bulk motion with respect to the 
CMB frame of 689 ± 178 km s~^ towards I = 343°, 6 = 52°. 
The 16 cluster sample studied here is much smaller than 
theirs, but the distance error per cluster is typically much 
smaller in this work. The results of comparing individual 
cluster distances for the 11 clusters in common are shown 
in the upper panels of Fig. |l^. The two sets of cluster dis- 
tances (peculiar velocities) disagree at the 99.9% (99.7%) 
confidence level, which indicates that the errors on one or 
both of the data sets are underestimated. Note that, in con- 
trast to the comparison between spiral and elliptical cluster 
samples, the mean redshifts and distances of our samples 
of cluster ellipticals, which typically include the BCG, are 
expected to be very similar to those of the BCGs alone (see 
the middle panels of Fig. p^ . 

The lower panels of Fig. |l^ show the comparison be- 
tween the clusters of this paper and BCG distances rederived 
by HE for the 9 clusters in common. The BCG distance indi- 
cator of HE is similar to that of LP but includes a correction 
to the BCG magnitude for the X-ray luminosity of the host 
cluster, because more X-ray luminous clusters tend to have 
brighter BCGs. The agreement between the our cluster dis- 
tances with those of HE is better than with those of LP, but 
is still marginal: the two samples disagree at the 94% con- 
fidence level. Part of the improved agreement results from 
the absence of A3381 from the X-ray sample, and part of the 
improvement is due to the X-ray correction itself: for the 9 
clusters in common, the reduced is 2.32 for the LP dis- 
tance indicator (still incompatible at the 98.7% confidence 
level) and 1.83 for the X-ray corrected BCG distance in- 
dicator. Nevertheless, there remain outlying clusters in the 
BCG samples: for both A262 and Coma (A1656) both the 
LP BCG distance indicator and the X-ray corrected BCG 
distance indicator of HE indicate peculiar velocities in excess 
of 2000 km s"^ (see Table ||). 
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Figure 11. Comparison between our cluster distances and peculiar velocities and literature values for the same clusters. Left-hand 
panels show the distance comparison. The dashed line has slope unity. The of the comparison and the number of degrees of freedom 
bottom right-hand corner. Middle panels show the comparison of cluster mean redshifts. Right-hand panels show comparison of cluster 
peculiar velocities. In all panels, the results of this paper are plotted on the horizontal axis. The upper panels show the results of HM, 
as rederived by Willick et al. (1997a); and lower panels show the results of Giovanelli et al. (1997a, b). 
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Figure 12. As in Fig. |ll| but note the expanded scale. Upper panels show the BCG results of Lauer & Postman (1994); lower panels 
show the X-ray corrected BCG distances of Hudson & Ebeling (1997). 
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Table 6. Peculiar velocity comparisons. Our measured cluster peculiar velocities axe compared with FP, TF and BCG results from the 
literature, for clusters in common. Literature sources are JFK96, HM data rederived by Willick et al. (1997a), Giovanelli et al. (1997b) 
(SCI), LP and HE. To the distance errors we have added in quadrature the error in the mean redshift, using the number of galaxies 
in the cluster and the cluster velocity dispersions derived in our work. For the BCG comparisons we adopt the mean redshift error of 
184 km quoted by Postman & Lauer (1995). Note that our results are very similar to those of JFK96, as expected since JFK is the 
principal data source for 6 clusters in this work. 

Fundamental Plane TuUy-Fisher Brightest Cluster Galaxy 



Hotel 


upp 


WJFK 










\^ KjlLLLLUSLHiO 




km 


km 


km 


km s"'^ 


km 


km 




Pisces 


131±208 




-136±229 


23± 343 








HMS0122 


-44±352 




104±298 


313± 291 








A262 


-259±340 




-783±355 


-34± 206 


-3119± 900 


-2953±1208 


see text 


A347 


-277±430 








-1467± 896 






A2199 


-342±325 






-308±1035 


464±1753 


40±1906 




A2634 


-960±364 




-844±590 


-156± 405 


437±1577 


-48±1428 


see text 


Coma 


0±204 


0±148 


-23±454 


231± 290 


2459±1372 


2051±1282 




A194 


743±230 


532d=188 






-1178±1088 


-1643± 914 




A539 


234±403 


213±285 


-410±604 




-225±1427 


-956±1376 




A3381 


578±593 


667±698 






-5905±2120 






A3574 


655±369 


556±388 




-131± 408 


20± 900 


711± 862 




DC2345-28 


-174±362 


-291±265 






-233±1467 


-1108±1481 




Hydra 


30±239 


87±261 


-85±322 


-453± 208 


-1063± 791 


-17± 584 





8 SUMMARY 

We have meeisured the mean peculiar motions of 103 early- 
type galaxies in 7 clusters in the PP region, and a further 
249 such galaxies in 9 calibrating clusters from the literature, 
using the inverse Fundamental Plane relation. This relation 
is found to have a distance error of 20% per galaxy. Our 
principal results are as follows: 

(i) Of 6 clusters in the PP ridge, none shows a significant 
motion with respect to the CMB frame. For the PP back- 
ground cluster J8, there is marginal evidence for 'backside 
infair into the PP supercluster. 

(ii) The PP supercluster has an insignificant net radial 
motion (— 60±220 km s~^) with respect to the CMB frame. 

(iii) An all-sky sample comprised of 16 clusters (with me- 
dian depth cz ~ 5500 km s^^) exhibits a bulk motion of 
420 ± 280 km s"^ towards I = 262°, b = -25°. 

(iv) Comparison of observed cluster velocities with pre- 
dictions from the IRAS 1.2Jy redshift survey yields f3i = 
OP'^ /bi = 0.95 ± 0.48, consistent with previous results. 

Our error analysis fully accounts for the uncertainties 
in the mean Hubble flow as well as the errors due to the 
merging of different spectroscopic systems. 

The bulk motion of the 16 cluster sample is consistent 
with the sample being at rest in the CMB frame, but is also 
consistent with the ~ 350 km s~^ motion found by Courteau 
et al. (1993). It is inconsistent with the ~ 700 km bulk 
motion found by LP. 

Our mean PP radial motion result is in apparent con- 
flict with the TF results of WUhck (1990, 1991) and HM 
which had found a ~ —400 km s^^ peculiar velocity of the 
PP ridge, but our results are in better agreement with the re- 
cent TF data of Giovanelli et al. (1996,1997a,b). Comparison 
of elliptical and spiral samples, is far from straightforward, 
however, due to the complex nature of the peculiar velocity 
field and the different regions probed by different surveys. 

The disagreement between our cluster distances and 
those of LP for the 11 clusters in common is statistically 



significant at the ^ 99.7% confidence (~ 3a) level indicat- 
ing that the errors of one or both of these data sets are 
underestimated. When the X-ray corrected BCG distances 
of HE are used, the disagreement is reduced to the ~ 94% 
(~ 2a) confidence level. 

The cluster peculiar velocity comparison is potentially 
the cleanest method of determining /? on linear scales be- 
cause it is relatively insusceptible to Malmquist and smooth- 
ing biases. However, a larger sample of clusters is required 
to reduce the random errors. 

The results from several large-scale peculiar velocity 
surveys are expected shortly. Of particular interest are the 
elliptical sample of the EFAR collaboration (Wegner et al. 
1996), the cluster TuUy-Fisher survey of Giovanelli, Haynes 
and collaborators (Dale et al. 1997) and the TuUy-Fisher 
field survey discussed by Strauss (1997; see also references 
therein for other peculiar velocity surveys). 

Finally, with collaborators, three of the present authors 
are currently obtaining FP data for ellipticals in an all-sky 
sample of approximately 50 clusters within a distance of 
12000 km s"^ This will allow a precise measurement of the 
bulk motion on these scales as well as a clean determination 
of /? in the linear regime. 



ACKNOWLEDGMENTS 

Michael Strauss is thanked for supplying to us the IRAS 
1.2Jy density field. This work has made use of Starlink facil- 
ities at Durham. JS and RJS acknowledge financial support 
from the PPARC. MJH acknowledges financial support from 
the PPARC; from a CITA National Fellowship; and from 
the Sciences and Engineering Research Council of Canada, 
through operating grants to F. D. A. Hartwick and C. J. 
Pritchet. 



REFERENCES 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy clusters in the Perseus-Pisces region — //. The peculiar velocity field 19 



Aaronson M., Bothun G., Mould J., Huchra J., Schommer R. A., 

Cornell M. E. 1986, ApJ, 302, 536 
Burstein D., Heiles C, 1982, AJ, 87, 1165 

Courteau S., Faber S. M., Dressier A., Willick J. A., 1993, ApJ, 
412, 51L 

da Costa L. N., Preudling W., Wegner G., Giovanelli R., Haynes 

M. P., Salzer J. J., 1996, ApJ, 468, 5L 
Dale D. A., Giovanelli R., Haynes M. R, Scodeggio M., Hardy E., 

Campusano L. E., 1997, AJ, in press 
Dckcl A., 1994, ARA&A, 32, 371 

Faber S. M., Burstein D., 1988, in Coyne G., Rubin V.C., eds. 
Proceedings of the Vatican Study Week, Large Scale Motions 
in the Universe. Princeton Univ. Press, Princeton, p. 135 

Faber S. M., Wegner G., Burstein D., Davies R. L., Dressier A., 
Lynden-Bell D., Terlevich R. J., 1989, ApJS, 69, 763 

Fisher K. B., Huchra J. P., Strauss M. A., Davis M., Yahil A., 
Schlegel, D., 1995a, ApJS, 100, 69 

Fisher K. B., Laliav O., Hofman Y., Lynden-Bell D., Zaroubi S., 
1995b, MNRAS, 272, 885 

Giovanelli R., Haynes M. P. 1985, AJ, 90, 2445 

GiovaneUi R., Haynes M. P., Myers S. T., Roth J. 1986, AJ, 92, 
250 

Giovanelli R., Haynes M. P., Wegner G., da Costa L. N., Freudling 

W., Salzer J. J., 1996, ApJ, 464, 99L 
Giovanelli R., Haynes M. P., Herter T., Vogt N. P., Wegner G., 

Salzer J. J., da Costa L. N., Freudling W., 1997a, AJ, 113, 22 
Giovanelh R., Haynes M. P., Herter T., Vogt N. P., da Costa 

L. N., Freudling W., Salzer J. J., Wegner G., 1997b, AJ, 113, 

53 

Gramann M., Bahcall N. A., Cen R., Gott J. R., 1995, ApJ, 441, 
449 

Guzman R., Lucey J. R., 1993, MNRAS, 263, L47 
Han M., Mould J. R., 1992, ApJ, 396, 453 

Huchra J. P., Geller M. J., Clemens C. M., Tokarz S. P., Michel 
A., 1993, the CfA redshift Catalogue, unpublished 

Hudson M. J., 1993, MNRAS, 265, 43 
Hudson M. J., 1994a, MNRAS, 266, 468 
Hudson M. J., 1994b, MNRAS, 266, 475 

Hudson M. J., Dekel A., Courteau S., Faber S. M., WiUick J. A., 

1995, MNRAS, 274, 305 
Hudson M. J., EbeUng H., 1997, ApJ, in press (HE) 
j0rgensen I., Franx M., Kjaergaard P., 1995a, MNRAS, 273, 1097 

{JFK95a) 

j0rgensen I., Franx M., Kjaergaard P., 1995b, MNRAS, 276, 1341 
(JKF95b) 

j0rgensen 1., Franx M., Kjaergaard P., 1996, MNRAS, 280, 167 
(JFK96) 

Lauer T. R., Postman M., 1994, ApJ, 425, 418 (LP) 
Lilje P. B., Yahil A., Jones B. J. T., 1986, ApJ, 307, 91 
Lucey J. R., Carter D., 1988, MNRAS, 235, 1177 (LC88) 
Lucey J. R., Gray P. M., Carter D., Terlevich R. J., 1991, MN- 
RAS, 248, 804 

Lucey J. R., Guzman R., Steel J., Carter D., 1997, MNRAS, 287, 
899 (LGSC) 

Lyndcn-BcU D., Faber S. M., Burstein D., Davies R. L., Dressier 

A., Terlevich R. J., Wegner G., 1988, ApJ, 326, 19 
Mathewson D. S., Ford V. L., Buchhorn M., 1992, ApJS, 81, 413 
Postman M., Lauer T. R., 1995, ApJ, 440, 28 
Raychaudhury S., 1989, Nat, 342, 251 

Santiago B. X., Strauss M. A., Lahav O., Davis M., Dressier A., 

Huchra J. P., 1995, ApJ, 446, 457 
Saunders W., et al., 1991, Nat, 349, 32 

Scodeggio M., Solanes J. M., Giovanelli R., Haynes M. P., 1995, 
ApJ, 444, 41 

Scodeggio M., Giovanelli R., Haynes M. P., 1997, AJ, 113, 101 
Shaya E. J., 1984, ApJ, 280, 470 

Smith R. J., Lucey J. R., Hudson M. J., Steel J., 1997, MNRAS, 
in press (PPI) 



Strauss M. A., Willick J. A., 1995, Phys Rep, 261, 271 
Strauss M. A., 1997, in Turok, N, ed., Critical Dialogues in Cos- 
mology, World Scientific, Singapore 
Tammann G. A., Sandage A., 1985, ApJ, 294, 81 
Watkins R., Feldman H. A., 1995, ApJ, 453, L73 
Wegner G., Haynes M. P., Giovanelli R., 1993, AJ, 105, 1251 
Wegner G., CoUess M., Baggley, G. Davies R. L., Bertschinger 
E., Burstein D., Mcmahan R. K. Jr., Saglia R. P. 1996, ApJS, 
106, 1 

Willick J. A., 1990, ApJ, 351, L45 

Willick J. A., 1991, PhD thesis. University of California, Berkeley 
Willick J. A., 1994, ApJS, 92, 1 

WiUick J. A., Courteau S., Faber S. M., Burstein D., Dekel A., 

1995, ApJ, 446, 12 
WiUick J. A., Courteau S., Faber S. M., Burstein D., Dekel A., 

Strauss M. A., 1997a, ApJS, 109, 333 
Willick J. A., Strauss M. A., Dekel A., Kolatt T., 1997b, ApJ, in 

press. 



© 0000 RAS, MNRAS 000, 000-000 



